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Abstract 

In this paper we discuss how the standard optimal Wiener filter theory can be 
applied, within a linear approximation, to the detection of an isotropic 
stochastic gravitational- wave background with two or more detectors. We apply 
then the method to the AURIGA-NAUTILUS pair of ultra low temperature bar 
detectors, near to operate in coincidence in Italy, obtaining an estimate for the 
sensitivity to the background spectral density of =10" 49 Hz -1 , that converts to an 
energy density per unit logarithmic frequency of ~ 8xl0~ 5 xp c , with 

p c «1.9xl0 - 26 kg/m3 the closure density of the Universe. We also show that by 

adding the VIRGO interferometric detector under construction in Italy to the 
array, and by properly re-orienting the detectors, one can reach a sensitivity of 
~ 6xl0" 5 xp c . We then calculate that the pair formed by VIRGO and one large 

mass spherical detector properly located in one of the nearby available sites in 
Italy can reach a sensitivity of ~2xl0 _5 xp c while a pair of such spherical 

detectors at the same sites of AURIGA and NAUTILUS can achieve sensitivities 
of «2xl0-6xp c . 



PACS 04.80.Nn, 95.55 Ym, 98.80.Es 



1. Introduction 

Resonant-mass gravitational-wave (GW) detectors are potentially 
interesting for the observation of an isotropic stochastic background of GW, as 
recently suggested by estimations 1 based on cosmological string models. These 
models predict for the dilatonic component of that background an almost 
frequency-independent GW spectral density Sh(co) that would put all the 
detectors under development on a comparable level of sensitivity. 

An experiment aimed at detecting a stochastic background of GW implies a 
correlation technique between the outputs of two or more detectors, possibly 
distant to each other to avoid spurious correlation of the detector noises. A sub 
optimal method has been discussed by Michelson 2 with reference to a generic 
detector, and a quasi-optimal method based on maximum likelihood estimation 
has been discussed 3 and applied to interferometric detectors. Also a linearized 
quasi-optimal method has been more recently discussed again in connection 
with an interferometer pair 4 . 

Here we first discuss the linearized quasi-optimal estimation method that 
can be obtained in a straightforward manner from the standard Wiener theory 
for optimal filtering. We show that it gives the same results as those of ref. 3 
for an array of detectors and, when applied to just a pair, of ref. 4, but in a 
more straightforward way. We then use the method to discuss the potential 
sensitivity of the pair of ultra low temperature resonant detectors AURIGA- 
NAUTILUS, which have been tested at two sites 5 at R ~ 400 Km apart in Italy. 

We then extend the calculation to the AURIGA-NAUTILUS-VIRGO 6 
array that will operate in Italy. 

We finally discuss a possible experimental detection strategy for spherical 
detectors and estimate the potential sensitivity of the pair formed by VIRGO 
and one spherical detector and of an array of two of such detectors. 

2. Isotropic GW stochastic background 

Isotropic stochastic background of GW has been widely discussed by many 
authors 7 . Here, we just summarise some of the results reported and give some 
useful details omitted in the literature. We assume that the metric tensor 
perturbation can be written as 




[1] 
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Under the assumption that the background is isotropic and stationary, the 
two-point two-time auto correlation of the metric tensor ^hy(r,t)h^(r' ,t')^ , 

where the brackets ( ) indicate an ensemble average, can be written as (see 
appendix A) 



(r , t)h kl (#•' , f )) = -L • J d0)S(0))e- iQ) ^- t ' } x 
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The functions 0Cq(x), OCi(x) and 0C2(x) in the above equation are given by 
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while the components Tjj k i of the matrices T n are given again in appendix A. 

12 
Note that, since Tmi=l, and r ^\\\\- r ^\\\\=^, and since 



lim cx Q (x) = 1, eqs. 2 and 3 give 
x— >o 



h 11 (r,t)h 11 (r,t')) = - 1 -- fdcoSCcOe-^') , 
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which is a standard normalisation in the context of signal analysis. 

The antenna pattern of any GW detector, in the approximation of detector 
size that is small compared with the gravitational wavelength, can be described 8 
by a symmetric traceless tensor D . By this we mean that the effective metric 

3 

perturbation sensed by the detector is h(t) = ^Dyhjj(t) with hij(t) the 

ij=l 

incoming wave. The explicit expression for the components Djj of D, for the 
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lowest longitudinal mode of a cylindrical GW antenna, is Dy = nm j - 8y / 3 , 

where n, with components ni, is the unit vector parallel to its axis. The lowest 
five degenerate quadrupole modes of a spherical detector have, instead, the 

/3 

following coupling tensors9 D™ (m=-2,..,+2) : D° =^-\2e+ z - e+ y ), 

D 1 =- X -e* z , D~ 1 = - X -e\, D 2 = ^e+ y and D~ 2 =-^e* y , where the 

tensors are defined in appendix A. 

Finally, for an interferometric antenna with its arms in the direction of the 
unit vectors n and m, we have Dy = (nmj - m^mj ) / 2. 

The cross correlation R^ b (t - t' ) = ^h a (t)h b (f between the effective 

metric perturbation signals h a (t) and h b (t) sensed by two detectors a and b 
located at r and r' , due to the stochastic GW background, can be written as 

oo 

Rj b (t - 1') = — ■ |dcoS(co)e" i(o(t " t ' ) x [ft oc (coR / c) + 
+ Clx ax (coR / c) + £l 2 ol 2 (coR / c)] 

3 

where Q. a = Dy D^T-j^ are three constants that only depend on the nature 
ijkl=l 

of the detectors, their locations, and their relative orientations. 

Let us now take two resonant bars. The force acting on antenna a placed in 



r a is 



, ,1 t d 2 h a (t) 
Fa(t) = -m a l a — [6] 

^ dt 



where m a is the effective mass (1/2 of the physical mass M) and l a is the 
effective length (for a cylinder l a =(4/7t 2 )L, with L the physical length of the 
cylinder). For two parallel oriented bars (the orientation that always maximises 
the correlation) we get 

(F a (t)iyt'))= 
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where: 



@(tf,x) = 



«p(x) 
8 



[7 + 4cos(2#) + 5cos(4#)] + 



[8] 



2 



[l + 2cos(2tf) + cos(4#)] + 



«2( x ) 
2 



[l-cos(4#)] 



here, $ is the angle between the direction of the detector axis (the same 
direction for both) and the straight line joining the sites where the detectors are 
located. The behaviour of the function x) is reported in figure 1 with $ 
as a parameter. 

It is worth pointing out that, if the detectors have both to be horizontal and 
still be exactly parallel, $ can only be $ = 90° . However, as in the following 
sections we will focus on observatories where R « R@, with R® the radius of 
the Earth, it is also worth noticing that in that case, whatever the value of 
the detectors can be made approximately parallel within an angle 8 = R/R® . 

For two spheres, each of the five degenerate quadrupole mode of the 
sphere acts as an independent detector 9 . We can write the GW force acting on 
the m-th mode as 



where now the effective mass and length refer to the mode (m a is the physical 

mass of the sphere and l a -0.3 of the diameter of the sphere 10 ) and h^(t) is 
obtained by contraction with the proper tensor D m . The cross correlation 
among the forces acting on the five modes of the two spheres separated by a 
distance R, assuming that all the modes are defined with respect to the same 
coordinate system, is then 




[9] 
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where 



S(m,x) 



-[a (x) + a 1 (x)] 

oc 2 (x) 
-[a (x)-a 1 (x)] 



m = 



m = ±2 



m = ±l 



[11] 



Note that for any value of m, E(m,x) has the same limit of 3/4 for and 
thus the sum over the five modes of two co-located spheres gives the usual 15/4 
factor coming from the solid-angle average of the energy emitted by a GW 
source. Figure 2 reports the behaviour of S(m,x) as a function of the 
dimensionless quantity x = coR/c. The correlation of the m = 0,±l modes 
decays faster than the m = ±2 ; therefore in a correlation experiment with two 
non co-located spheres (coR / c > 5) we have only these two modes relevant for 
the detection of the GW stochastic background. 

3. The linearized estimation method 

The stochastic force resulting from the GW background adds just an extra 
contribution to the Gaussian noise of a single detector. As a consequence, a 
measurement with a single detector has to rely on an a priori estimate of the 
detector noise that has to be subtracted from the data, a procedure whose 
accuracy is always highly questionable. 

A much safer method is obviously to cross correlate the outputs of two or 
more detectors with no common source of noise. Correlation of noise sources, 
such as electromagnetic interference, seismic noise, etc., is expected to decay 
with the distance between the detectors. Cross correlation between detectors 
located far apart is then the most advisable procedure to estimate the 
background. 

The Wiener theory of optimal estimation can be applied both to 
deterministic signals and to stochastic processes buried in detector noise. In this 
section we show that through a suitable small-signal linearization, a Wiener-like 
method can be applied to a cross-correlation estimate of the amplitude of a 
stochastic background signal of known power spectrum driving a set of noisy 
detectors. 
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Assume then that the power spectrum of the GW background S(co) in 

eq. 2 can be written as S(co) = A S(co), where S(co) is a known function of 

co, while the amplitude A is unknown and has to be estimated. S(co) can be a 
flat function of CO, as suggested by some recent string theory 1 , or the power 

law S(co) CO ^, which is predicted by the standard inflationary cosmology 7 . 
The output data streams of a set of N detectors can then be written: 

x a (t) = r| a (t) + As a (t) 1 < a < N , [12] 

where r| a (t) are stochastic processes that describe the noise in each detector 

and s a (t) is the signal due to the stochastic GW background. 

The statistical properties of the Gaussian processes in eq. 2 can be 
summarised as follows: 

(r[ a (t)) = (s a (t)) = 

(r| a (t)ri b (t')) = R a (t-t')5 ab l<a,b,<N 

[13] 

(s a (t)s b (f)) = R ab (t-f) 
(n a (t)s b (f)} = 

The correlation R ab (t - t' ) among the GW signals available at the output 

S 

of the detectors are 

oo oo 

R ab (t - f ) = Jdt"H a (t") J"dt'"H b (t w )Rf (t - t' -t" + t w ) [14] 



where R^ b (t - t' ) is related to S(co) by the same eq. 5 that relates R ab (t - t' ) 
to S(co). The response functions H a (t) of the a th detector translate the input 
signal h a (t) into the output signal 

oo 

As a (t) = |H a (t / )h a (t-t')dt / . [15] 
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We now look for an optimal estimator of the spectrum amplitude by 
writing down the most general bilinear combination of the outputs of N 
detectors: 

N j j 

A2= Ll T dt I T dt 'gab(t,f)x a (t)x b (f) , [16] 

a,b=l 

where the "filter" functions g ab (t,t') have to be chosen such that A is an 
unbiased estimator: 

A 2 ^ = A 2 [17] 

and that its variance o^ 2 =^A 2 j ^-^A 2 ^> is minimal. 

These requests lead in general to a non-linear problem. However, we show 
in appendix B that if 

As a (t)«n a (t), [18] 

which is likely to be the case for actual detectors, then the problem can be 
linearized, and g a b (t,f ) obeys the integral equation 



[19a] 



T T 

2 J dt" J dt'"g ab (t", r )R* (t - t" ) R b (f - i'") = - 1 R s ab (t - f ) 
-T -T 

for a^b; while 

g ab (t-t') = fora = b [19b] 

The Lagrange multiplier X/2 in eq. 19a is obtained from 
N 



T T 

J dtj dt'g ab (t,t')Rf(t / -t) = l 



[20] 



a,b=l 

and we can also show that: 
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al 2 =-- [21] 

A 2 2 

Eq. 19a can be solved numerically. However much of the information 
about the solution can be obtained by assuming that g a b(t,t') decays rapidly as a 
function of lt-t'l. To be more specific, if the optimal bandwidths of the detectors 
are of order 1/t and if the light travel time is ~R/c, it is then reasonable to 
assume that g a b (t, t' ) — > if {|t - 1' | — R / c}/x » 1 and then choose to 
integrate on a time 2T»2t+R/c. If this is the case, we can take T« °° in eq 19a 
and solve it in terms of Fourier transforms: 

to («>)—- S ° b((0 K ■ [22] 
where 



oo 

gab (t, f ) « g ab - f ) = - 1 - f g ab (coy^dCO, 
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[23] 



— oo 



and where we have defined (co) and S^ b (co) as the Fourier transforms of 

II S 

R* (t - t' ) and R* b (t - t' ), respectively. 

To evaluate the variance of the estimate, we can then substitute this 
solution in eq. 20: 



-2 
a A2 = 



+oo +oo T T 

S s ab (co)sf(co') 



= Y^ f fdcodco' 8 ^ 8 ^ [ fdtdt'e^'H—'^ 
4^(2ti) 2 J J 2S*(co)S b (co) J J 



a*b v ' -oo-oo nv / nv / _ T _ T [2 4] 

+oo +oo 

=y^r fdcodco- 5 - w t (co ' )sm[(co " CQ,)T] 5 T (co-co-) 

^2tc J J S*fco)S b fco) (co-co')T 



S£(co)S°(G)) (co-co')T 

d^D — oo — oo 

T 

<? , x 1 f ioot i 1 sin(Tco) . _ . . . . 

where o j (co) = — I e dt = is a fmite-time approximation to the 

2tt J 71 co 

-T 

Dirac 8 function, which reduces to 5(co) in the limit T — > oo . If we assume 
again that the observation time T is large enough to have null correlation 
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functions for |t-t'|>T, then 8j(co) is a sharply peaked function in a very 



small frequency range compared to the scales on which the functions S ab (co) 

S 

and S a (co) change. In this case then eq. 24 gives 



°A 2 = 



+ 00 



LIS 

a<b -00 



dco 



Sf(co) 



SS(co)sS(co) 



-1/2 



[25] 



In eq 25 we have also changed from the sum over a^b to that over a<b, 
which contains half the terms. Note then that, if we define 



+00 



t r 

71 J 
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Sf(co) 



-1-1/2 



SS(co)sE(co) 



[26] 



which is the uncertainty of the estimate that only uses the data from detectors a 
and b, eq. 25 becomes 



-=T- 



^2 

a A2 



a<b 



^2 
a ab 



[27] 



The result in eq. 25 can be further clarified: let us call H a (00) the transfer 
function of the a th detector, i.e. the Fourier transform of the response function 

H a (t) appearing in eq. 15. Then clearly 



S ab (co) = S aD (co) H a (co) H D (co) , where S aD (co) is the Fourier transform 



,ab 



ab 



of R^ b (t-t') in eq 15. By defining the equivalent input detector noise as 
S b (co) = (co) / H a (co) 2 , eq. 25 can be rewritten as: 



^2 

°A 2 = 



LIS 

a<b -oc 



dco 



Sf(co) 



Sj(a>)s8(a>) 



-1 



[28] 



or, by using eq. 2 and the definition of S(co) 
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G A 2 = 



OO 

LIS 

a<b -o< 



dco X 



X 



{s(co)[^ oc (coR / c) + Q.i «! (coR / c) + £l 2 a 2 (<° R / c )]} 

S^(co)sE(co) 



-i-l 



[29] 



Note that eq. 27 agrees with the results of ref. 3, which have been derived with 
a different approach. 

The weight function g ab (t - t' ), which is just the optimal bilinear filter 
for the search of the stochastic background, can be explicitly evaluated by 
solving eq. 19a in the Fourier space: 

gab(a>) = <4 2 x 

S(co)[£2 oc (coR / c) + £2j a x (coR / c) + Q. 2 oc 2 (coR / c)] P°] 

S^co)S^(co) 

The filter depends on the shape of power spectrum S(co), but the dependence 
should be very weak for any detector array that includes resonant detectors. In 

fact such detectors have comparatively narrow bands, so if S^(co) is its 

equivalent noise input spectrum, then S^(co) grows very rapidly outside the 
post detection band centred at some centre frequency COq in the kHz range. As a 
consequence, the template in eq. 30 decreases rapidly outside the same band. In 
this limit, S(co) can be approximated by a constant S(co) ~ S(c0q )• 

If, instead, the array includes only wide bandwidth detectors, a set of 

filters with S(co) °c co~^ should be constructed and the data processed with the 
different choices of y predicted by different cosmological models. The filters 

9 

that maximise the signal-to-noise ratio A / o ~ 2 gives then an estimate of y. 
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4. Implications for AURIGA, NAUTILUS, VIRGO and large mass 
spherical detectors. 



In this section, we use eq. 29 to discuss the sensitivity of the AURIGA and 
NAUTILUS pair (soon to be operated in coordinate coincidence) to the GW. 
stochastic background. We also discuss the expected performance of the 
AURIGA-NAUTILUS -VIRGO array and the sensitivity of arrays including one 
or more large mass spheres. 

Let us consider the simplest model for a resonant antenna, in close analogy 
with the viewpoint suggested by Giffard 11 . The a th antenna is considered as a 
simple harmonic oscillator of mass m a and length l a with resonant angular 

frequency co a and decay time T a , excited by the force. The position x a (t) of 
the oscillator mass is read by a suitable position transducer. The transfer 
function from the input force to the output displacement is 

1 co 2 

H a (co) = ^ 5 y [31] 

2m a co a - co +ico/x a 

Within this model the oscillator is driven into random motion by the sum of the 
brownian noise force and the back-action noise force of the position transducer. 

The position transducer also contributes with an additive position noise x a (t). 

Both the total noise force f n (t) and the additive displacement noise x a (t) are 

assumed to have white spectra with values Sf and S a respectively. Sf and S| 
can be parametrized as: 

Sf = N a hk a 

N a h [32] 



S a = 
°x - 



k a 
K n 



where N a represents the noise quantum number, which has to be N a > 1/2 , and 



kn = A/S? / i s a n °i se "stiffness". These quantities have an immediate physical 
meaning because they directly relate to the antenna burst sensitivity h^ and 
post detection bandwidth Aco^ respectively. Here, h^n is defined as the 
amplitude of the burst, of centre frequency =co a and duration =l/co a , which 
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gives a signal to noise ratio of one. The relation of N a and to h^ n and 
Aco a d is: 



'mm 



1 



2N a h 



1/2 



v m a co a j 



k a 

; Aco a pd ~^ 



m a co a 



[33] 



The total noise at the antenna output is given by the sum of a narrow band and a 
wide band contributions: 



m 



a(cO^-C0 2 ) +(C0/T a ) : 



[34] 



This noise can be referred to the antenna input, as if it were a spectrum, 
dividing S a (co) by the squared antenna transfer function. The resulting noise 
spectrum is shaped as the inverse of a lorentzian function and can be written as 



S a h (o>) = 



Aco 



P d 



CO 



[35] 



where cof = cof 

a a 



a 2 4 



~ co a is the frequency at which the noise reaches 



a minimum and which coincides in practice with the resonance frequency of the 
antenna. 

SKco) reaches its minimum value at CO ~ co a : 



[sg(co)l . «SE(© a )« 

L Jmin 



(2h a nin )^ 



co: 



Aco pd - 



8k B T 



m a la»aQa 



[36] 



where Q a =co a T a and T* is the temperature of the antenna. 
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4.1 Two bars 



Let us consider the case of two bars. Using the physical length L a and mass M a 
eq. 36 gives for a single bar: 



Sh(co a )«1.3xlO" 45 x 



X 



At V 2 f 



3m 



2300 Kg 



50 mK 



lKHz 



V V a J 



uo 7 ; 



Hz 



-1 



[37] 



To evaluate the sensitivity of a pair of bars to the stochastic background we 
substitute eq. 35 in eq. 28, obtaining: 



I 



X 



Aco pd Aco^ d (4h min h min )^ 



co 8 S 2 (co)© z (#,coR/c)dco 



X p ^ 2 ~1 T 2 

(co^-co 2 ) +(coAco; d ) (co 2 -co 2 ) +(coAcoj; d ) 



[38] 



To carry out our calculations, we now assume that the detectors are 
parallel and that they have the same sensitivity, the same post detection 
bandwidth but, to take into account the real physical situation, we allow for 

slightly different resonant frequencies. For (cOa/AcOpdj »1, the integral in 
eq. 38 can be approximated by 



C A 2 



[«1 



— CO i 



2TAco; d 



© 2 (i3,coR/c) 



1/2 



[39] 



where co = (co a + C0b) / 2 is the mean frequency of the two antennae and where 

we have taken S(co = co a ) = 1. 

For R = 400 Km, which is the distance of the AURIGA-NAUTILUS pair, 
and with an average frequency of co = (2n • 920) rad/ s, coR/c = 7.7 while the 

expected sensitivity is \sl] ~ 1.7 X 10" 45 Hz" 1 . Taking # = 90° one gets: 

L J min 
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C^2 « 10 



-49 



( 2T > 


-1/2 




/Aco pd ^ 


v lyear , 






30 

y 



-1/2 



Hz" 



[40] 



provided that the detuning of the resonance frequencies is small compared with 
the optimal bandwidth. This last condition is, for instance, matched with 25 % 

approximation if |co a - <x> b | /AcOpd < 1 / 2. With an effective bandpass > 20 Hz, 

this implies that the two detectors have to be matched within 10 Hz, which 
appears to be feasible. 

However, in order to have both detectors parallel to each other and as 
parallel as possible to the other cryogenic detectors already in operation 12 , 
AURIGA and NAUTILUS are presently oriented with r> = 52°. This value 
gives a sensitivity of roughly a factor of 2 worse than that in eq. 40 (Fig.l). 

To reorient the detectors is technically feasible, but it is doubtful whether 
the factor of 2 is worth the loss of parallelism with the remaining detectors. 



4.2 Two bars and one interferometer. 

We now discuss the potential sensitivity of the AURIGA-NAUTILUS - 
VIRGO array, assuming for VIRGO the planned position (Lat = 10° 30' E, 
Long = 43° 40' N) and orientation (one arm at 26°, the other at 296°). There 
are many noise sources in an interferometric antenna that have been carefully 
estimated 13 . However, the resulting total noise power spectrum in the frequency 
range of interest here, i.e., around 1 kHz, is dominated by the shot-noise 
contribution: 



s; irgo (co)~s virgo 



( \ 2 

CO 



ho 



v^oy 



400 Hz < CO / 27i< 4000 Hz, [41] 



where, from the published curves of the expected noise spectrum 13 , one can 

estimate S^ rgo = 1.6 x 10" 45 Hz" 1 if co is taken co = 27t920 rad / s. This gives 

a noise value close to the minimum value of the input noise of AURIGA or 
NAUTILUS. 

One can apply eq. 5 to the pairs AURIGA- VIRGO and NAUTILUS- 
VIRGO, which have intersite distances of ~ 220 Km and ~ 260 Km 
respectively, corresponding to co a R/c = 4.2 and co a R/c = 5 at a frequency 



15 



of 920 Hz. The correlation value corresponding to the AURIGA- VIRGO pair 
turns to be «2 greater than that for NAUTILUS -VIRGO. 

The sensitivity to the stochastic background for each bar-interferometer 
pair and for one year of integration is obtained by substituting eqs. 40 and 35 in 
eq. 28 and by integrating: 



( Toa 1 ovirgo 



o k2 = 



[s a h ] . s 

L Jmin 



virgo 
ho 



TAc0p d @ 2 (coR/c) 



[42] 



where we can see that the overall bandwidth is set by that of the resonant 
antenna. The resulting sensitivity for the pair AURIGA- VIRGO is 

C =2xlO- 49 Hz" 1 , assuming for AURIGA co a /Au)p d = 30 and the present 

orientation; while for the NAUTILUS VIRGO pair, again with the present 

orientation, 0-2 = 3.5 x 10~ 49 Hz _1 . For the array as a whole, using eq. 27, 

these values give 0^2 = 1.3 x 10~ 49 Hz -1 . With the bar detectors oriented with 

$ = 90°, a choice that maximise the overall sensitivity, one gets instead 

a 22 = 8xl0" 50 Hz _1 . 

Although the improvement in sensitivity with respect to the pair of bars 
is almost negligible, such a three-detectors detection would be of paramount 
importance in ruling out spurious effects. 



4.3 Two spheres 

Finally we study the performances of two spherical detectors, which have 
been recently proposed as possible next-generation resonant antennas. 

Eq. 38 only involves the detectors post detection bandwidth and input 
noise. We can easily compare the sensitivity of a bar to that of a sphere made 
of the same material, and with the same resonant frequency, as 14 : 



o bar 
5 h 



= 1.17 
Jmin 



q sphere 
5 h 



J mm 



(M s /M b ), [43] 



where M s and M^, are the physical masses of the sphere and the bar 
respectively. For a sphere of diameter D, made of the same material of the 
present bars (Al 5056) we find 
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Sh(co a )«7xlO""' x 



x 



3m 



38000 Kg 



50 mK 



( 



lKHz 



\3 



V V a J 



uo 7 ; 



Hz" 



[44] 



Under the same approximation used for eq.38, the sensitivity of a sphere pair, 
for each mode of the sphere, can be written as: 



O k 2 



^Lnl^Ln^d^-^bl/AcOpd)^ 



nl/2 



2TAC0pd 



E 2 (m,coR / c) 



[45] 



If we consider two co-located spheres, the overall sensitivity is higher by a 
factor V5 with respect to that for a single mode, since we can add the signals on 
the five modes of the sphere. If the spheres are far apart, fig. 2 shows that only 
for the m = ±2 modes is the decay of the correlation with distance as slow as 
that for the bars. For the other modes, the correlation decays much faster, so 

for x>5, the overall sensitivity is higher only by a factor of ~ a/2 with respect 
to that for a single mode. 

The expected sensitivity of two co-located 3-m-diameter sphere, made of 
Al 5056, is 



G^ 2 « 4 x 10 



-52 



lyear 



Qy A copd 

30 



v-l/2 



Hz 



-1 



[46] 



if one takes the weighted average of all the modes. As all the parameters in eq. 

9 

45 are the same for the five modes except for the coefficients S (m,coR / c), 
the error on the weighted average is given by the same eq. 45 by replacing 

2 5 

S (m,coR/ c) with the sum £E(m,coR/c) which is reported in Fig. 3. For 

m=l 

two spheres at the AURIGA and NAUTILUS sites, coR/c = 7.7, one would get 



then o 1 2 -2x10 

A z 



-51 



2T 
lyear 



-1/2 



30 



,-1/2 



Hz 



-1 



This figure may improve by changing the sphere material and/or 
increasing the sphere diameter. For instance if one is able to fabricate two 4-m- 
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diameter copper alloy sphere (250 tons), the above reported sensitivity reaches 
about 4 10-52 Hz- 1. 



4.4 One sphere and one interferometer 

In fig. 4 we show the correlation function for one interferometer and a 
sphere. The function is obtained by summing up the contribution coming from 
all the sphere modes. If the sphere and the interferometer are not at the same 
site, then the function depends on the angle 9 between one of the interferometer 
arm and the line joining the two sites. For a sphere and an interferometer like 
VIRGO located at the same site, the figure gives: 



o~2 « 8 x 10 



-51 



2T 
lyear 



-1/2/ 



co a /Aco pd 



n-1/2 



30 



Uz 



[47] 



As the site of VIRGO is fixed, one can try the exercise to locate a sphere in 
one of the three major laboratories available close to VIRGO. Those are the 
AURIGA and NAUTILUS sites and the large underground laboratory of Gran 
Sasso 15 . For these laboratories coR/c and 9 take the values 
coR/c = 4.2, 9 = 4° ,coR/c = 5, 9 = 68.7° and coR/c = 5.3, 9 = 90° 

respectively. With T=l year and COa/AcOp^ = 30, these figures give sensitivities 

that are respectively : o^ 2 ~ 2.5 x 10" 50 Hz _1 for the AURIGA site, 

a~ 2 ^xlO" 50 !^" 1 for that of NAUTILUS and a~ 2 =6xl0" 50 Hz _1 for 
the Gran Sasso Laboratory. 

5. Conclusions 



The sensitivities reported in this paper can be expressed in term of the 
ratio ^gw(co) of the mass-energy density per unit logarithmic frequency of the 
GW stochastic background to the closure density p c of the universe. In fact the 
spectrum of the stochastic background S (co) can be written as 16 : 

. 16G71 2 _ . . 1A _46^ , Y2x Jixl kHz^ 3 
S(co) = 3— £2 GW (co)p c =8xlO 4D £2 GW (co) 



or 



co j 



An experiment involving at least one resonant detector, bar or sphere, 
operating at 1 kHz would than measure ^gw(co) with an uncertainty: 
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^ GW ~8xl(T 6 



G A 2 



10 



-50 



where p c = 1.9xl0 -26 Kg/m 3 has been assumed. The sensitivities discussed in 
sec. 4 can be then recast in terms of ^gv(co) according to table I. 





Table I 




Estimated sensitivities of various detectors arrays. 


Detectors array 


°A 2 ( HZ ) 


"GW 


AURIGA-NAUTILUS 


2x10-49 


2x10-4 


present orientation 






AURIGA-NAUTILUS 


10-49 


8x10-5 


best orientation 






AURIGA-NAUTILUS- 


1.3x10-49 


10-4 


VIRGO 






present orientation 






AURIGA-NAUTILUS- 


8x10-50 


6x10-5 


VIRGO 






best orientation 






VIRGO and one 38 ton 


2.5x10-50 


2x10-5 


sphere 






Two 38 tons spheres at 


2x10-51 


2x10-6 


AURIGA and 






NAUTILUS sites 







These figures have to be compared with experimental limits already 
existing. Microwave background measurements by COBE 17 limit £2gw to about 
10" 8 but just at very low frequencies (<10"15 Hz). Pulsar timing observation also 
give limits < 10 -9 but again for very low frequencies ( ~ 10 -8 Hz)* 8 . An 
estimate from earth quadrupole oscillation 19 in the mHz range gives £^GW< 3. 

Direct experimental measurements on ground in the kHz range have been 
performed in the past with sensitivities significantly lower than those indicated 
in table I. A pioneering measurement with split bars 20 gave a limit of £2gw< 

3xl0 3 in the kHz range. More recent estimates from the background noise of a 
cryogenic detector 2 ! gi ve £2q W < 3x10 2 . 

Table I then shows that already in the near future experiments involving 
resonant mass detectors can provide unprecedented upper limits to the GW 
stochastic background at kHz frequencies. Future sensitivities of bars- 
interferometer and spheres experiments may hope to go near or beyond the 
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limit put by nucleosinthesys considerations at £2gw~10~ 5 and to detect such 
fundamental cosmic signals. 
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Appendix A 



The coefficient of the expansion in eq. 1 can be written as a superposition of in- 
going and out-going monochromatic plane waves, 



h ij (A:,t) = hi] > (A:)e- lkct +h^(A:)e- 



[Al] 



where we have denoted by k the modulus of k . Since the metric has to be a 
real function, we have hjj(fc,t) = hjj(-fc,t), which implies hy (k) = hy (-k). 

We then assume that the amplitudes \i^(k) and h-j~(A:) are Gaussian 

stochastic processes with zero mean and a given power spectrum and that the 
integral in eq. 1 exists 22 . 

As each plane wave is transverse and traceless, in the coordinate system 

(x,y,z) where k = k^sinScos^M^ + sin sin §u y + cos9« z ), we can write 



hjj*(*) = hjj* (*) + hjj* (*) = 

= A ik (9,^h^ + (A:)eJ m +h^ X (A:)eJ m jA m ;(e,rt 



where ejj - 



(\ 0^ 
0-10 
0y 



and 



(o i o^ 

1 

y 



[A2] 



are the two independent 



polarisation tensors in the TT frame and Ay(9,(|)) is the Euler matrix that 

rotates the axes, from the TT system, where k = (0,0, k ), to the (x,y,z) frame. 
It follows that 



hij (*) = 



= h^ "(ft) 



+ 



2 2 2 I 2 \ 

cos 9 cos ()) — sin <j) I cos 9 + 1 1 sin <j) cos § 

(2 \ 2 2 2 

cos 9 + 1 1 sin <j) cos <|) cos 9 sin <j) - cos <|) 

- sin 9 cos 9 cos <|) - sin 9 cos 9 sin <j) 



sin 9 cos 9 cos § 
- sin 9 cos 9 sin <j) 
sin 2 9 



+ 
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X 



^-cos9sin2(|) cos9cos2(|) 
cos9cos2(|) cos9sin2(|) 
sin 9 sin (() - sin 9 cos <|) 



sin 9 sin <j) ^ 
- sin 9 cos <]) 




[A3] 



_|_ ^ 

i.e., hjj > (A:) = h _> (A:)\|/J (9,(|)) + (£)\|/ x (9, <|>) with obvious definitions of 
the v^'s. 

In order to satisfy the reality condition for the transform in eq. 1 we need: 

h^~(A:) = h^ + (-£)\|/jj*(7t-9,4) + 7i) + h^ X (-k)\\ffj(n-Q,^ + ll), or, as 
can be easily checked from eq. A3 by substituting 9 — » 71 - 9 and <|) — » 7t - <|) , 



x* 1 

h-j~(A:) = h _> (-£)\|/J(9,(|)) - h - ^ (-k)\\f x (9, <\>), so that only two complex 

I w 

amplitudes are really independent and we call them simply h (A;) and h (k). 

Assume now that the processes h + (k) and h x (k) are uncorrected, 
stationary, and isotropic: 



h + (Jfc)h x V)) = lh + (k)h x *(k')) = 0, 



h + (A:)h + *(A:')) = (h x (A:)h x *(A:')) = (27r) 3 S(k)5(k-k')5 2 (^, MA: ') 



[A4] 



where S(k) = S(IA;l), and are the unit vectors parallel to k and k' and 

8 2 is the Dirac 8 function on the unit sphere. This implies that the 
two-point, two-time correlation function of the GW stochastic background 
(|hjj(r,t)h lm (r' ,t')^ is just a function of the modulus of the distance between 
the two points and of the modulus of the time difference 



hii(r,t)h lm (r\t')) = 



x 



x 



HI cos[kc(t - f )]{[v5(*)vL(*) + Vij (^)v t /r m (^)]s(k)e 1 ^- r ')}d 3 k = 



[A5] 
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= ^-3 • |dkk 2 S(k)e" ikc(t " t ' ) x J* d(f> J* sin 9d9 x 
^ 2n ^ — 

x ¥ ;(9,« ¥ f m (9,« + ^ 

where we have assumed r - r' = R^sin 9' cos (j)' u x + sin 9' sin u y + cos 9' u z j . 

8 k 2 

By taking (|) -0, 9=0 and S(co) = S(k) , the angular part of the integral 

1571 c 

in eq. A5 can be performed explicitly and eq. 2 is obtained with the following 
values for the matrices TO, T 1 and T 2: 



T = 



' e + 
K xz 

2 xy 


f e + - 
K xz 



±e x 
2 *y 



+ 



+ 



'xz *xy 







e + -2e + 
K xy ^ K xz j 



'Jty 



±e x 
2 xy 





±e x 
2 xy 

K xz 




'jcy 






-2e 



+ 
xz 



2 












> x 
yz 



x A 

xz 

X 

yz 

, 





f° 





n 




(0 







e x = 
K xz 











e x = 

' c yz 








1 




vl 











1 


o, 



and is the null 3x3 matrix. 



[A6] 







(I 





0^ 




ri 





0> 




(0 


1 




where 


K xy 





-1 





e + = 











' K xy 


1 












K.0 









K.0 









^0 





o, 
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Appendix B. 

By substituting eq. 16 in eq. 17 we immediately get 
N 



£| dt|dt'g aa (t,t')R*(t-t') + 



a=l 

N 

+ 

a,b=l 



T f T 
dt d 
T J-T 

T f T 
dt d 
T J-T 



JN j T 

Y,[ dt|dt'g ab (t,f)R s ab (t-f) = 



[Bl] 



from which eqs. 19b and 20 readily follow. 

The ^(a 2 j ^ term which enters the estimation of o^ 2 can be calculated by 
using 19b as 

T rT rT rT 



\( A2 ) ) = L J T dt J T dt/ J T dt// Jdrgab(t,t')gcd(t // ,0 

a^b,c^d 

x(x a (t)x b (t')x c (t")x d (r) 



X 

T J-T J-T J-T " TO01 

a*b,c*d . [B2J 



Then using eq. 13 and the known rule for zero mean Gaussian random variables 

(xyzw) = (xy)(zw) + (xz)(yw) + (xw)(yz) we get 

(( a 2 f ) - ( a 2 f = £ |_ t t dt -£ T T dt"|df" gab (t, t') gC d(t", n ■ 

\K(t" - t)R b (r - t')5 ac 5 bd + K (r - t)R b (r - t')5 ad 5 bc + 
+A 2 [Rf(r-t)R b (r-t')5 bd +R ad (r-t)R b (r-t')5 bc + [B3] 
+ R bd (r - t-)Ri(t" - t)5 ac + R bc (r - t-)Ri(r - t)5 ad ] + 

+A 4 [Rf(t" - t)R bd (t w - f) + R ad (t w - t)R bc (t" - t')]} 

Assume now that the correlation signal As a (t) is negligible in comparison to 

the intrinsic noise T| a (t), as is the case in the real physical situation. Hence, in 
eq. B3 we can neglect the terms containing A 2 and A 4 , getting 
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°h =(( kl f ) - ( kl f - 2 L DlVlVlV' x 



a*b T J " T - 1 M 

Xgab(t,t')gab(t",t'")R a n (t"-t)Rl;(t'"-t-) 

The problem reduces to a constrained variational problem where, with the help 

9 

of the standard lagrangian multiplier technique, we minimise 0^2 under the 
constraint of eq. Bl. The functional 

•T rT pT pT 

dt dt' dt" dl 
-X J-T J-T J-T 



Mg ah ,k) = Y\2\ T dt [ T dt'[ T drfdrx 

™~ J— T J— T J— T J— T 

a^b L 

x[gab(t,t')gab(t",r)Rn(t"-t)R^(r-t')]+ [B5] 

+x\ dt f dt' gab (t,t')R s ab (t'-t) 



reaches its minimum when 

= l<a,b<Nanda*b, [B6] 



8A(g ab ,X) 



5 Eab 

i.e., when eq. 19a is obeyed 
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Caption to Figures 



Fig 1 Correlation function @(9,x) for two parallel bar detectors as a function of 
the reduced distance x=coR/c, with R the detector distance. The different curves 
are parametrized by the value of the angle 9 between the detector axis and the 
line joining the detector sites. The flattest curve corresponds to 9=0, the curve 
of maximum oscillation is for 9=71/2. The dotted line corresponds to 9=1 rad, 
which is the orientation of the AURIGA -NAUTILUS pair. 

Fig 2 Correlation function E(m,x) for the -2< m <2 modes of two spheres as a 
function of x=coR/c, with R the detectors distance. The different curves refer to 
the different modes according to the legend in the insert. The modes relate to a 
reference frame where the z axis is along the line joining the detectors 



interferometer and a sphere as a function of the reduced distance x = coR/c . 
Here the function <&(m,x) is the correlation function between the 
interferometer and the m th mode of the sphere. 



Fig. 3 The overall correlation function 




function of the reduced distance x = coR/c. 



Fig. 4. The total correlation function A(x)= ( m > x ) for 




an 
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Fig. 1 
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Fig. 2 
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x= coR/c 



Fig. 3 
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Fig. 4 
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